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ABSTRACT 

Gravitational lensing induces significant errors in the measured distances to high- 
redshift standard candles and standard sirens such as type-la supernovae, gamma-ray 
bursts, and merging supermassive black hole binaries. There will therefore be a signif- 
icant benefit from correcting for the lensing error by using independent and accurate 
estimates of the lensing magnification. Here we investigate how accurately the magni- 
fication can be inferred from convergence maps reconstructed from galaxy shear and 
flexion data. We employ ray-tracing through the Millennium Simulation to simulate 
lensing observations in large fields, and perform a weak-lcnsing reconstruction on the 
simulated fields. We identify optimal ways to filter the reconstructed convergence maps 
and to convert them to magnification maps, and analyse the resulting relation between 
the estimated and true magnification for sources at redshifts zg = 1 to 5. We find that 
a deep shear survey with 100 galaxies/arcmin 2 can help to reduce the lensing-induced 
distance errors for standard candles/sirens at redshifts zs ~ 1-5 (z§ rj 5) on average 
by 20% (10%), whereas a futuristic survey with shear and flexion estimates from 500 
galaxies/arcmin 2 yields much larger reductions of 50% (35%). For redshifts z$ ^ 3, a 
further improvement by ~ 5% can be achieved, if the individual redshifts of the galax- 
ies are used in the reconstruction. Moreover, the reconstruction allows one to identify 
regions for which the convergence is low, and in which an error reduction by up to 
75% can be achieved. Such strongly reduced magnification uncertainties will greatly 
improve the value of high-redshift standard candles/sirens as cosmological probes. 

Key words: gravitational lensing - distance scale - gravitational waves - super 
novae: general - gamma-ray burst: general 



> 1 INTRODUCTION 



5-H 



j^ As suggested by Lemaitre (19271 and confirmed by Hub- 



ble (19291, there is a correlation between the distances and 



redshifts of galaxies that arises from the expansion of our 
Universe: galaxies at greater apparent distance show more 
redshifted spectral features. The exact relationship between 
distance and redshift changes for different cosmological theo- 
ries or different cosmological parameters. Accurate measure- 
ments of the distance-redshift relation may therefore help to 
discriminate between models and to constrain cosmological 
parameters ( Linder| 2008). 

Distance estimates based on the apparent magnitudes 
of type-la supernova (SN) have provided constraints for the 
current expansion rate and the mean matter density of our 
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Universe, and substantiated the evidence for the currently 
favoured ACDM model (e.g. |Komatsu et al.|2009[ ), since SN 
observations indicated a cosmological constant or another 
dark-energy component is currently accelerating the cosmic 
expansion ( Riess et al.|1998 Perlmutter et al.|1999 1 . In con- 
junction with other cosmological probes, current SN data 
even puts constraints on the properties of the dark energy 
( Kowalski et al.||2008| [Rubin et al!]|2009 l. 



While type-la SN can only be seen out to redshifts z ~ 
2, gamma-ray bursts (GRB) can be detected at much higher 
redshifts. There is growing evidence for correlations between 
observed and intrinsic burst properties (e.g. |Ghirlanda et al.| 
2006 1, with which GRBs can provide distance estimates out 
to redshifts z > 6 ( Schaefer|2007 |. When added to distance 
data at lower redshift, the higher-redshift bursts help to 
break parameter degeneracies in cosmological models and 
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make GRBs a valuable discriminant between dark-energy 
models ( |Schaefer|20"07||Liang et al.|2008||Amati et al.|2008 |. 
Another promising way to obtain very accurate distance 
measurements in the future are observations of gravitational 
waves emitted by merging supermassive black hole binaries 



pact objects ( Rauch 1991 Metcalf and Silk 1999 Mortsell 



(MSBs) ( |Schutz|[l986| |Holz and H ughes 2005| ). From the 
gravitational wave phasing, it is possible to determine the 
luminosity distance to the source and the redshifted mass, 
M(l + z), to high precision. However, the redshift cannot 
be measured from the gravitational waves alone, so the use 
of MSBs as distance probes will rely on the identification 
of an electromagnetic counterpart to measure the redshift 
of the merger event (Lang and Hughes 2008). Detection of 
the gravitational-wave signals from MSBs is one of the main 
goals of the planned Laser Interferometer Space Antennae 
(LISA). MSBs are very luminous and can be observed out 
to the cosmic horizon. Current estimates suggest that LISA 
will see a few tens of events per year, of which most will 
have distance measurements to less than ~ 10%, and of the 
order of two per year will have distance measurements to 
less than 1% and sky localisation to within a square degree 
( Arunet al.|2009b| >. 

The distance measurements for standard candles like 
type-la SNe, GRBs, and standard sirens like MSBs are based 
on the comparison of the electromagnetic or gravitational ra- 
diation energy emitted by the source (assumed to be known 
at least statistically) and the radiation flux or amplitude 
received by the observer's detector. Apart from the emit- 
ted signal strength and the distance between source and 
observer (and other factors), the received signal strength 
is also affected by gravitational lensing: Gravitational de- 
flection by intervening matter structures may focus or de- 
focus the radiation towards the observer, resulting in an 
increased or decreased detector signal. This gravitational 
(de-)magnification induces an additional scatter in the dis- 
tance estimates (see, e.g., Holz and Linder|2005 and refer- 
ences therein). 

Evidence for gravitational (de-)magnification has been 
found in the GOODS SN d ata (| Jonsson et al.|20"06| ) and the 
Supernova Legacy Survey ( Jonsson et al.|2010b l. The lens- 
ing effect is small compared to other sources of uncertainty 
for current SN samples (e.g. |Riess et al.]|2004[ ), but it will 
become significant for large future high-redshift (i.e. with 
z > 1.5) SN samples (Holz and Lindcr 2005 Jonsson et al 



|2010b[ ). In particular, pencil-beam surveys might be compro 
mised by gravitational lensing because of the large spatial 
correlation of the magnification ( Cooray et al.| [2006). 

There is also strong evidence that high-redshift GRBs 
are significantly affected by gravitational lensing and magni- 
fication bias (e.g. |Porciani et al.||2007| |Wyithe et al.||2010[ |, 
which may lead to large biases in their deduced distances 
if not properly accounted for. Moreover, magnification scat- 
ter will almost certainly be the dominant source of error in 
distance measurements using MSBs ( Arun et al.|2 009a). 

While a source of noise for distance measurements, mag- 
nification of standard candles/sirens has also been consid- 
ered a useful signal of the cosmic matter distribution. The 
observed scatter in SNe brightnesses can be used to con- 
strain the amount of dark matter in the form of massive com- 



1 http://lisa.nasa.gov/ 



et al.|200T Minty et al. 2002 I, or to measure the power spec- 
trum of density fluct uations on small scales ||Metcalf| 1999 \ 
and on larger scales (Dodelson and Vallinotto 2006). Cor- 
relations between galaxies and SN brigtnesses can be used 



to constrain the properties of galaxy halos (Metcalf 2001 
Jonsson et al. 2010a|c| 



When measuring distances to standard candles/sirens, 
one can try to correct for the gravitational lensing effects 
by inferring the magnification from other observations. For 



example, Gunnarsson et al. ( 2006 1 and Jonsson et al. ( 2009 1 



suggested modelling the matter distribution between the ob- 
server and a source using a model for the halos of the ob- 
served galaxies near the line of sight. These foreground mod- 
els could reduce magnification scatter for high-redshift SN 
by up to a factor of 3. However, the success of this technique 
relies heavily on the validity of many assumptions about the 
detailed relation between observed light and total matter. 



Dalai et al. ( 2003 1 used convergence power spectra to es- 



timate the correlation between the true convergence (which 
is closely linked to the magnification) and the convergence as 
reconstructed from background galaxy shear measurements. 
They found that even with a very high number density of 
galaxies for the convergence reconstruction, the reduction in 
the magnification-induced scatter was small. Based on sim- 
ilar calculations, |Shapiro et al.| ( |2010[ | estimated that con- 
vergence maps reconstructed using galaxy flexion in addition 
to shear might help to reduce the lensing-induced distance 
errors by up to 50%, making this particular 'delensing' ap- 
proach promising. 

In this work, we use ray-tracing through the Millen- 



nium Simulation (Springel et al. 2005) to investigate how 



accurately the lensing magnification of standard candles and 
standard sirens in a ACDM universe can be inferred from the 
shear and flexion measured in the images of distant galaxies. 
The accuracy with which the magnification can be estimated 
directly determines the lensing-induced uncertainty in the 
distances of type-la SNe, GRBs, and MSBs that remains 
after delensing, i.e. after a correction for the magnification 
effects has been applied. In a companion paper (Gair, King, 
& Hilbert, in prep.), we will discuss the consequences of 
delensing for the use of gravitational-wave sources as cos- 
mological probes in comparison to other probes. 

Compared to the magnification correction scheme sug- 



gested, e.g., by |Gunnarsson et al. (20061, the method dis- 
cussed here has the advantage that it does not require any 
assumptions about the relation between the distribution of 
visible and dark matter. This method consists of (i) recon- 
structing the convergence toward a SN, GRB, or MSB as 
accurately as possible from the shear and flexion measure- 
ments of galaxy images near the line-of-sight, (ii) estimating 
the magnification from the reconstructed convergence, and 
(iii) correcting the observed SN, GRB, or MSB signal using 
the inferred magnification. 

As will be shown below, for the correction to yield sat- 
isfactory results, one must reconstruct the convergence on 
sub-arc-minute scales (see also Dalai et al.|2003 \ , where the 
non-Gaussian nature of the underlying density field becomes 
important. Moreover, a realistic description of the noise of 
the convergence reconstruction on these scales has to take 
into account the discreteness and randomness of the galaxy 
image positions. Finally, the relation between the magni- 
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fication and the noisy reconstructed convergence is non- 
trivial. Such complications usually have to be neglected in a 
treatment relying on convergence correlations and Gaussian 
statistics (as was done, e.g. by Dalai ct al. 2003 Shapiro 
et al. 20101. In contrast, all these effects can be included 



in the ray-tracing simulations used here. This approach not 
only provides more realistic predictions about the perfor- 
mance of magnification correction schemes, but also an op- 
timised magnification estimate, and more detailed informa- 
tion about the residual magnification after correction. In 
particular, the simulations yield detailed probability distri- 
butions for the magnification and its residual as a function of 
the reconstructed convergence. Such information is valuable 
as input for Bayesian parameter estimation for cosmological 
models, and helps to further reduce lensing-induced errors 
and biases on cosmological parameters ( Hirata et al.|[2010 
Shang and Haiman|2010| 



The paper is organised as follows: In Section [2] we dis- 
cuss gravitational lensing and how it affects the distance 
estimates of standard candles and standard sirens. In Sec- 
tion [3] we describe the methods for simulating noise-free 
lensing maps and for creating shear and flexion maps with 
realistic noise properties. Moreover, we discuss how the in- 
formation in the noisy maps is used to estimate the local 
magnification. Results for the magnification correction are 
presented in SectionH] The paper concludes with a summary 
and discussion in Section [5] 



2 THEORY 

2.1 Gravitational lensing 

The photons and gravitons (and any other particle species, 
e.g., neutrinos) from a distant source may be deflected gravi- 
tationally by intervening matter inhomogeneities before they 
reach the observer. The deflection, called gravitational lens- 
ing (see, e.g., Schneider et al.|2006 for an overview), causes 
shifts in the observed image position relative to the 'true' 
source position, i.e. the sky position at which the source 
would be seen in the absence of gravitational lensing. Here, 
the lensing effects are described by the lens mapping, 



£ 



[0, 



(0, 



(3 = f3(0,z), 



(1) 



which relates the angular image position — (6-1,62) of a 
point-like source at redshift z to its true (unlensed) angular 
source position /3 = (/?i,/?2)Fl 

The distortion induced in the images of extended 
sources like galaxies can be quantified to first order by the 
Jacobian of the lens mapping: 



A(0,*) = 



(dpi(0,z)\ 
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(2) 



This equation defines the convergence k(0,z), the shear 
7(0, z) = 71 + i72, and the asymmetry cj(0, z). Higher-order 
image distortions can be quantified by higher derivatives 

2 The flat-sky approximation is used here to avoid irrelevant 
complications due to the spherical geometry of the 'real' sky, and 
is valid given the survey areas that we consider. 



of the lens mapping, which can be combined, e.g., into the 
spin-1 flexion 

(3) 
(4) 

(5) 

The (unsigned) image magnification is given by: 
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2.2 Lensing effects on distance estimates 



(6) 



The lensing magnification influences the observed signal 
strength of distant SN, GRBs, and MSBs, and thereby the 



inferred distances (e.g. Gunn|[l967 Schneider and Wagoner 



19871 jKantowski et al.||1995J |Frieman||1996| |Wambsganss 



et al. 19971. In a nutshell, the radiation power Pg emitted 



by a standard candle/siren at redshift zs and the power Pn 
received by a detector are related byfj 



c nf i(0,z s ) 

"r — -7^0 — 1 — r " s - 

AumOs) 



(7) 



Here, cr denotes a detector-specific constant, fi(0,zs) the 
lensing magnification in the direction of the image, and 
Aum the luminosity distance. Hence, 



A U m(zs) = \/cbju(0, ZS 






(8) 



An error Sfj, in the magnification causes an error SDi u 
in the luminosity distance: 

5-Dlum ~ — 7: -Dlum, 

2 (A 

or in logarithmic form, 



(9) 



<51n(Aum) ~ --<51n(/i). 



(10) 



2.3 Quantifying the magnification uncertainty 

Delensing methods provide estimates fi cs t of the true mag- 
nification /1 of a standard candle/siren. The resulting error 
of the delensing method is encoded in the distribution of the 
residual magnification /i rcs = fi//j, cs t. A worthwhile delens- 
ing method should result in a residual magnification whose 
statistical dispersion is substantially smaller than the dis- 
persion in the true magnification fi. 

In order to quantify the performance of delensing meth- 
ods, one needs a way to quantify the dispersion in the resid- 
ual magnification /i rcs and compare it to the dispersion in 

3 In practice, received radiation energies or photon numbers oc 
Pr are measured for type-la SN and GRBs, while gravitational- 
wave detectors record wave amplitudes <x \/Pr.. This difference 
does, however, not affect the resulting relation between inferred 
distances and magnification. 
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the uncorrected magnification fi. A common way to quan- 
tify the dispersion of a distribution is the standard deviation. 
However, for point-like sources, the magnification distribu- 
tion Pfj,(n) oc /i -3 for fi ^> 1 and so the standard deviation of 
the magnification is formally divergent. In this work we will 
therefore quantify the dispersion of the magnification using 
the distribution of the logarithm of the magnification 



A = ln(/x) 



(11) 



and its standard deviation, <j\, which is finite for point-like 
sources. Similarly, the dispersion of the residual magnifica- 
tion /^res will be quantified by the standard deviation <J\ tee 
of its logarithm 

A rcs = ln(jUres). (12) 



2.4 Magnification estimates from galaxy shear 
and flexion 

There are various methods to obtain an estimate for the 
magnification (J,(0, z) for a given line of sight and redshift 
z. One can, for example, try to infer the matter distribution 
near that line of sight from the observed distribution of mass 
tracers like galaxies, and then calculate the expected mag- 
nification fi(0, z) from the inferred matter distribution (e.g. 



Gunnarsson et al. 2006 1. One can also try to infer the magni- 



fication from source number density counts (e.g. Broadhurst 
et al.p995 Zhang and Pen|2005 1 , or the image sizes of ex- 
tended sources (e.g. Bartelmann and Narayan|1995[ ). 

Here, we focus on methods that use weakdensing recon- 



struction ( Kaiser and Squires 1993 Bacon et al. 2006 1 to 



obtain magnification estimates (see also Dalai et al. 2003 
|Shapiro et al.|[2010[ ). These methods 'reconstruct' the con- 
vergence from the shear and/or flexion measured in the im- 
ages of distant galaxies and use the result to calculate the 
magnification. In their conventional form, they neglect the 
asymmetry, and assume that the convergence, shear, and 
flexion are well approximated by derivatives of a common 
lens potential ip{6, z)y\ 

(13a) 
(13b) 

(13c) 

(13d) 

(13e) 

Hence, the convergence can be reconstructed from the shear 



4 These approximations are exact for pure E-mode lensing. B- 
mode lensing can be taken into account by taking a second po- 
tential into consideration (e.g. |Schneider et all]|2002| |Hirata and| 
Scljak 2003). B-modes are, however, usually very small compared 
to E-modes ||Jain et al.|2000||Hilbert et al.|2009i. 
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in Fourier space. Here, hats denote Fourier transforms, i.e. 
1 



k(l,z) = —l d A cxp(-ie-0)K(0,z), 

£ — {£\,£\) denotes the two-dimensional Fourier vector, I — 
(-1 + i^2 denotes the complex Fourier wave number, and an 
asterisk denotes complex conjugationrl The reconstructed 
convergence is then translated into a magnification estimate 
via an approximate version of Eq. Jm applicable in the limit 
of weak lensing: 

Mo st(0, 2 )~ [1-k(0,z)]~ 2 (15a) 

&1 + 2k(0,z). (15b) 

The magnification estimate (J,est(0, z) is then used to 'delens' 
the standard candle/siren. 

There are several complications in practice. First, nei- 
ther the shear 7, nor the flexion T and Q can be directly 
measured, but only their reduced versions = 7/ (1 — «), 



F = X7gg, and G = Vg g ( Schneider and Seitz|1995 Schnei- 
der and Er||2008 l. This complication can be taken into ac- 
count, e.g., by iterative reconstruction methods (Seitz and 
|Schneider||1995p . Moreover, 7 w g, T » F, and Q w G for 
most parts of the sky. It is thus assumed in the following 
that galaxy shapes provide estimates for the shear and flex- 
ion (and we postpone the proper treatment of reduced shear 
and flexion to forthcoming studies). 

Even from a perfect shear or flexion map, the conver- 
gence can only be reconstructed up to a constant because 



of the mass-sheet degeneracy (Falco et al. 19851. In addi- 



tion, the boundary conditions influence the reconstruction 
near the edges of finite fields. These problems can, however, 
be alleviated by using a large field of view for the recon- 
struction (whose mean convergence is close enough to zero, 
as discussed by |Shapiro et al.||2010[) and special fin ite-field 
inversion techniques (e.g. Seitz and Schneider 1996). 

Furthermore, the galaxies that can be used to measure 
the shear and flexion have a finite number density n ga i and 
thus only provide shear and flexion estimates at a finite num- 
ber of discrete and randomly distributed points. Moreover, 
the shape of a galaxy image provides a very noisy estimate 
for the shear and flexion. Therefore, spatial interpolation 
and averaging are necessary to obtain a convergence esti- 
mate with acceptable signal-to-noise in the desired direction 
of a standard candle/siren. The smoothing erases real struc- 
ture on small scales, which causes additional deviations of 
the reconstructed convergence from the 'true' convergence. 
The optimal smoothing scale is determined by a trade-off 
between the aim to reconstruct as much small-scale struc- 
ture of the convergence as possible and to average over as 
many galaxy images as possible, so as to reduce the noise 

5 The definition of the Fourier transform used here differs from 
the one used, e.g., by Kaiser and Squires 1 1993 1. 
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due to the dispersion in their intrinsic shapes (Dalai et al 



2003) 



Finally, the galaxies do not all lie at the redshift zs 
of the standard candle/siren, but have a certain redshift 
distribution p ga i(2 ga i) depending on their intrinsic redshift 
distribution and on the depth of the survey. A reconstruction 
scheme not using individual galaxy redshiftsj only yields a 
noisy and smoothed version of the effective convergence: 



ft e ff(0) 



dz ga lPgal(2gal) k(0, %al) 



(16) 



If the redshift distribution p g ai(%ai) is not sharply peaked 
around the redshift zs of the standard candle/siren, the ef- 
fective convergence K e g(0) may deviate substantially from 
the convergence k(8,zs), which may cause a systematic as 
well as a statistical error in the reconstruction of k(6,zs). 
Spatial variations of the 'local' redshift distribution of the 
galaxies constitute an additional source of noise. 



3 METHODS 

The complications discussed above limit the accuracy with 
which the magnification can be estimated from the weak- 
lensing reconstruction. Therefore, accurate estimates of their 
effects are needed in order to evaluate the performance of 
the weak-lensing reconstruction. Since it is difficult to obtain 
accurate estimates for all these effects by analytical calcula- 
tions, we use a numerical approach. 

Using ray-tracing through structure-formation simula- 
tions, we create a suite of simulated fields of view populated 
with galaxy images at random positions with noisy estimates 
of the local shear and flexion. The simulated galaxy images 
are then used to reconstruct the convergence in the simu- 
lated fields, and the reconstructed convergence is used to 
calculate an estimate of the magnification. The 'true' mag- 
nification (obtained from the ray-tracing) in the field is then 
corrected by the magnification estimate from the reconstruc- 
tion, and the resulting residual magnification is analysed. 



3.1 The ray-tracing 



Here, the Millennium Simulation (MS) by Springel et al. 
1 2005), a large high-resolution iV-body simulation of cosmic 
structure formation, is used to simulate the matter distri- 
bution between observer and source. The cosmological pa- 
rameters for the MS are: a matter density S1m = 0.25 and a 
cosmological-constant energy density £Ia = 0.75 (in units of 
the critical density), a Hubble constant h — 0.73 (in units 
of 100kms -1 Mpc -1 ), a primordial spectral index n — 1 
and a normalisation parameter as = 0.9 for the linear den- 
sity power spectrum. The simulation was run using a par- 



allel TreePM version of GADGET2 (Springel 2005) with 
iV p = 2160 3 particles of mass m p = 8*6 x 10 s h M© in a 
cubic box with comoving side length L — 500/i -1 Mpc, and 
an effective force softening length e = 5/i -1 kpc. 



The gravitational lensing effects from the matter struc- 
tures in the MS are calculated by the Multiple-Lens-Plane 
ray-tracing algorithm described in Hilbert et al. ( |2009j ). The 
algorithm takes into account the lensing effects of the dark- 
matter structures on scales > 5/i -1 kpc, which are rep- 
resented by the simulation particles in the MS. In addi- 
tion, the effects of the stars in galaxies with stel lar masses 
^ 10 9 /i -1 Mq (as predicted by the model of De Lucia 



Hilbert et al. 



and Blaizot (2007l) are taken into account as described in 



(2008) 



The multi-plane algorithm is used to simulate 32 fields 
of view, each with an area of 4 deg x 4 deg (providing 512 deg 2 
total area) and covered with a regular grid of 4096 2 pixels 
(providing a resolution of 3.5arcsec) in the image plane. For 
each pixel, a ray is traced back through the MS by the al- 
gorithm to calculate the magnification, convergence, shear, 
and flexion for sources on each of the 44 lens planes that 
span the redshift range from z — to z = 5.7. 

Every pixel of the ray-tracing fields (excluding a 
15 arcmin margin, which might be affected by artefacts in 
the weak-lensing reconstruction, and strong-lensing regions 
with det(A) ^ or tr(A) ^ 0) is used as a sample line-of- 
sight towards a standard candles/siren. Here, we consider 
only volume-limited samples of standard candles/sirens. 
Thus, each pixel is weighted in the statistical analysis by 
its inverse magnification to account for the magnification 
bias ( |Canizares|1981| |Hilbert et al.|2007b l. 



3.2 The simulated weak-lensing surveys 

Simulated catalogues of lensed galaxy images with a given 
number density n ga i are generated with random sky posi- 
tions uniformly distributed in the fields of view. The red- 
shifts of the simulated galaxies are drawn from a distribution 



that matches observations (Brainerd et al. 1996) 



PgaH%alj 



P 



Zgjd 
20 



exp 



%alY 

20 ) 



(17) 



where a, 0, and 2o are parameters (with zq proportional to 
the median redshift of the simulated survey). 

Here, we consider three idealised weak-lensing galaxy 
surveys: an advanced survey, a futuristic survey, and a per- 
fect survey. For the advanced survey, we assume a galaxy 
density n ga i = 100 arcmin -2 and a redshift distribution with 
a = 2, P = 3/2, and median redshift z me dian = 1-5 appropri- 



ate for a deep space-based lensing survey (cf., e.g., |Schrab- 
back et al.||2010"| . For the futuristic survey, we assume a 
very high density n ga i — 500 arcmin -2 of usable galaxy im- 
ages having a redshift distribution with a = 0.8, — 2, 
and 2 me dj an = 1.8 appropriate for an ultra-deep survey ( |Coe 
et al.|[2006 | with extremely high spatial resolutionlj In ad- 
dition, we consider a 'perfect' survey, whose noise- free con- 
vergence maps are directly created from the ray-tracing. 

For the advanced and futuristic surveys, we consider 
both the case that only shear estimates are used, and the 



6 Taking into account the individual galaxy redshifts (e.g. in a 
tomographic reconstruction) complicates the reconstruction, but 
may improve the magnification estimates. A simple way of using 
the individual redshifts to improve the reconstruction is discussed 
in Section|4.5| 



7 Even with future telescopes, such high densities of galax- 
ies with shape measurements might only be reached for a very 
deep pointed observation, which could then be supplemented by 
a wider survey with lower number densities I Shapiro et al. 12010 1 . 



© 0000 RAS, MNRAS 000, 000-000 



6 S. Hilbert, J.R. Gair, & L.J. King 



case that shear estimates are combined with flexion esti- 
mates. The shear and flexion estimates of the simulated 
galaxies are obtained from interpolating the shear and flex- 
ion from the ray-tracing onto the galaxy positions and red- 
shifts. The shape noise in each galaxy image is simulated by 
adding a random number drawn from a normal distribution 
with standard deviation a 1 al = 0.2, a? a] = 0.5arcmin , 
and ag al = 0.9arcmin _1 , to each component of 7, J-, and 
Q respectively (Rowe et al., unpublished, see Shapiro et al 



20101 



3.3 The weak-lensing reconstruction 

For the weak-lensing reconstruction, we consider a mesh- 
based method that allows one to employ Fast Fourier Trans- 
forms (FFT). Each field of view is covered with a regular 
mesh of 4096 2 pixels (i.e. the same mesh geometry as for 
the ray-tracing). In each pixel, the shear 7 is estimated by 
weighted sums over the shear estimates 7^ of the simulated 
galaxy images in the field: 



«„,« 



i gal 'gal 
E* ^gal 



.» 



The weights w ' Y are determined 



*B 



(i) 
"gal 



Mo-o%) 



(18) 



(19) 



where 6 denotes the sky position of the pixel centre, the 0^ 
denote the sky positions of the galaxy images, 



wg(6) = exp 



261 



(20) 



and # w = 1.75arcsec denotes a weighting scale, which is 
chosen slightly smaller than the pixel size. The shear es- 
timate on the regular mesh is then Fourier transformed 
by FFT methods ( |Frigo and Johnson||2005 1 , and converted 
into a convergence estimate k 7 using Eq. ( 14a I and setting 
k^(£ = 0) = 0. If only shear information is used, k 1 is used 
directly as an estimate of the 'raw' convergence in Fourier 
space prior to additional smoothing: 



fiJraw(,£j — K^\K). 



(21) 



In the case that flexion information is used as well, es- 
timates kjr and kg of the Fourier-space convergence from 
the galaxy flexion T and Q are computed in a procedure 
analogous to that for k 1 . The convergence estimates from 
shear and flexion are then combined into a joint estimate by 
( Okuraet aL|2007| ): 



,m 



^ + m 2 °rL + w 2 °gL 



(22) 



The smoothing inherent in the averaging ( 18 \ may be 



insufficient to reach an acceptable signal-to-noise ratio of 
the convergence estimate in real space. The raw convergence 



estimate K ra w is thus smoothed by a Gaussian low-pass filter 
with filter scale 8 S , i.e. 



Acst (£) = exp 



/"'■r. 



Al), 



(23) 



before the resulting convergence estimate is transformed 
back to real space by FFT methods (j 

3.4 Magnification estimates from noisy 
convergence maps 

The convergence map generated from the weak-lensing re- 
construction can be converted into a map of estimated 
magnification in several ways. If the estimated convergence 
ftcst(#) is close to true convergence k(0, zs), and if the shear 
and rotation are small, 



A cst (0) = -21n[l-K cst (6>)] 



(24) 



provides a good estimate for the logarithmic magnification 
X(0, zs) of a standard candle/siren at position and red- 
shift zs Fj This simple estimate might however fail, if the 
estimated convergence K cst (0) differs substantially from the 
convergence n(G, zs), e.g. because of shape noise or because 
the weak-lensing survey probes a redshift range very dif- 
ferent from the redshift zs of the standard candle/siren. In 
certain cases, the simple estimate may perform even worse 
than simply ignoring lensing and assuming A cst = (see 



Section 4.2 1 



If the conditional distribution P\\ KBBt (A|K es t) of the true 
magnification A for a given estimated convergence «: cs t is 
known, one can construct a bias-free magnification estimate 
that minimises the dispersion <7A rcB in the residual magnifi- 
cation by (see Appendix [A[) : 

AeetW = <A) AKst [««t(0)], (25) 

where (A) A i K (re C st) is the expectation value of the true 
magnification A for a given estimated convergence K cst : 

(26) 



See Section 4.5 for using redshift-dependent weights. 



< A >A| Kcat (««»0 = / dAp A | Koat (A|K cs t)A. 



This estimate can be computed from the lensing simulations 
by spatial averages over the simulated fields. The estimate 
is optimal in the sense that no other magnification estimate 
based on the estimated convergence yields a smaller disper- 
sion in the residual magnification. For example, the resulting 
dispersion in the residual magnification is never larger than 
the dispersion in the uncorrected magnification (i.e. assum- 
ing A est = 0). 



4 RESULTS 

4.1 Perfect reconstruction 

We first consider the case that one can somehow obtain a 
perfect estimate of the true convergence k(0, zs) for sources 



See Section |4.4| for using Wiener filters instead of Gaussian 
smoothing. 

10 One might consider using the full Eq. Kh instead, but tests 
show that this does not improve the magnification estimate for 
realistic noise levels in the reconstruction (except near strong- 
lensing regions where the discussed reconstruction method is not 
suitable anyway). 
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Figure 3. The dispersion 0V, of the residual A res (dashed line) 
for sources at redshift 2g = 1.5 as a function of the filter scale 8 S 
used to generate the convergence estimate re es t by convolving the 
convergence re(2g) with a Gaussian filter. The solid horizontal line 
marks the magnification dispersion A res = A without correction. 



Figure 1. The joint distribution of the logarithmic magnification 
A = ln(^t) and the convergence re for sources at redshift zg = 
1.5 (lighter areas correspond to higher probability densities on a 
logarithmic scale). 
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Figure 2. The probability distribution PA ros (A r es) of the residual 
logarithmic magnification A rcs = A — A cs t for sources at redshift 
2g = 1.5. The residual distribution for the optimal estimate J25J 
with the convergence to the same source redshift (dashed line) is 
compared to the case without correction (i.e. A GS t = 0, solid line). 



at the same position and redshift Zs as the standard can- 
dle/siren. This case provides hard limits to what can be 
gained by delensing in more realistic cases. 

The joint distribution p\^ K {(j,, k) of the logarithmic mag- 
nification A = ln(/i) and convergence k, for sources at redshift 
zs — 1.5 is shown in Fig. [T] The joint distribution is almost 
degenerate with the probability concentrated around the re- 
lation A = — 21n(l — k). This indicates that one can get a 



very good magnification estimate from the Eq. (25 1, if the 



convergence can be reconstructed with sufficient accuracy. 
Moreover, the optimal estimate (p5l is virtually identical to 



the simple estimate ( 24 1 



timate is illustrated in Fig. [2] The distribution of the un- 
corrected residual A rcs = A for standard candles/sirens at 
redshift zs = 1.5 has a standard deviation crA rcs = 0.075. 
The corrected residual A rcs = A — A cst using the optimal 
estimate (251 and the convergence n cst (9) — k(0,Zs) has a 



The effect of delensing using a perfect convergence es- 



fifteen times smaller standard deviation o"a ios = 0.005. The 
small but non- vanishing residual dispersion is due to the fact 
that the magnification also depends on the shear, but this 
dependence is ignored in the estimate |25|Fj 

Even with perfect measurements for individual galaxies, 
a magnification estimate will make use of galaxies spread 
over a certain patch of the sky. This amounts to averaging 
or smoothing the convergence over the area occupied by the 
galaxies. The impact of smoothing on the dispersion of the 
residual magnification is shown Fig. [3] The magnification 
correction is already substantially degraded (reaching 50% 
of the uncorrected dispersion), if a smoothed version of the 
convergence with a filter scale S — 10 arcsec is used. This 
illustrates the need to faithfully reconstruct the convergence 
on scales of a few arc seconds, if one wants to decrease the 
dispersion in the residual magnification by a factor two com- 
pared to the uncorrected magnification. 



4.2 Realistic reconstruction 

As a more realistic scenario, we assume that the conver- 
gence K e st is reconstructed from the shear measured in an 
advanced survey with a galaxy density n ga i = 100 arcmin -2 
and median redshift z mc di an = 1-5. The joint distribution 
of the logarithmic magnification A for sources at redshift 
zs = 1.5 and the estimated convergence K ost is shown in 
Fig. p] The joint distribution roughly follows the relation 
Acst = — 21n(l — West), but there is a large scatter around 
it. Moreover, the mean (A),, (n BB t) as a function of the 
estimated convergence K es t deviates substantially from that 



11 Tests show that taking the shear dependence into account 
does not improve the magnification estimate noticeably in the 
presence of realistic noise. We thus refrain from including an ex- 
plicit shear dependence. 
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Figure 4. The joint distribution PA,re cat (A, ftest) of the loga- 
rithmic magnification A for sources at redshift 23 = 1.5 and 
the convergence K e st reconstructed from the shear of galaxy im- 
ages in a survey with galaxy density n ga i = lOOarcmin - 2 , me- 
dian redshift z me dian = 1-5, and a Gaussian smoothing with 
filter scale 9 B = 20 arcsec (lighter areas correspond to higher 
probability densities on a logarithmic scale). Also shown are 
the relation A = — 21n(l — re BS t) (solid line) and the relation 
■* = Wx\K, cst ( K ost) (dashed line). 



relation, in particular for k, < 0. This shows that the es- 



timate ( 24 1 may differ substantially from the optimal esti- 
mate ( |25[ ) for a realistic noise level. 

Fig. [5] illustrates how much the dispersion <7A rca hi the 
residual magnification A rcs of a standard candle/siren at 
zs — 1.5 can be reduced with an advanced survey for partic- 
ular choices of the filter scale 9 S . For 9 S ^ 30 arcsec, both es- 



timates ( 24 1 and ( 25 1 yield a similar reduction in the residual 



dispersion. For 9 S < 30 arcsec, however, the optimal estimate 
( 25 1 performs significantly better than the simple estimate 
(24f. 

The optimal estimate always yields a residual disper- 
sion 0"a ics smaller than the uncorrected dispersion a\, but 
the simple estimate results in o~x raB > cr\ for 9 S < 10 arcsec 
if only shear data is used. In that case, the noise in the 
reconstructed convergence is much larger than the signal 
(i.e. the true convergence). The simple estimate translates 
this noise into a large noise in the estimated magnification, 
which thus acquires a dispersion exceeding the dispersion of 
the true magnification. In contrast, the optimal magnifica- 
tion estimate responds more weakly to the noisy convergence 
and thus keeps the residual dispersion below the uncorrected 
dispersion. 

If only shear data is available, the residual dispersion 
(TA r eB f° r sources at zs = 1.5 is reduced to 81% of the uncor- 
rected dispersion o\, when the simple estimate and an opti- 
mal filter scale 9 S w 30 arcsec is used. The optimal estimate 
yields a slightly smaller residual dispersion tJ\ It , s — 0.79<ta at 
the optimal filter scale 8 S w 25 arcsec. Further improvement 
is obtained if flexion information is available in addition to 
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Figure 5. The dispersion <TA roa of the residual logarithmic mag- 
nification A rcs = A — A cs t of sources at redshift 23 = 1.5 as a 
function of the filter scale 8 B used to reconstruct the convergence 
Kest from an advanced survey with n ga i = lOOarcmin -2 . Re- 
sults are shown for the simple estimate l |24[ ) (dashed lines) and 
for the optimal estimate ]25[ (dotted lines) with convergence re- 
constructed from shear (a) and from shear and flexion (b). The 
solid horizontal line indicates the dispersion A rcs = A without 
correction. 



the shear. In this case, the residual dispersion is reduced to 
76% by the simple estimate and a filter scale 9 S m 15 arcsec. 
The optimal estimate ( 25 1 performs best at 9 S sa 5 arcsec, i.e. 



at the mean separation of the galaxy images. The residual 
dispersion is then 70% of the uncorrected dispersion. 

In Fig.[6l the probability distribution of the uncorrected 
magnification is compared to the distribution of the residual 



magnification using the optimal estimate ( 25 1 for correction 



The residual distribution is not only narrower, but also less 
skewed than the distribution of the uncorrected magnifica- 
tion, making it a better candidate for a Gaussian approxi- 
mation. 

In a statistical analysis, one might prefer to use the indi- 
vidual probability distributions PA roa |fc OBt (Ares, i\K BS t,i) of the 
residuals X rcs ,i appropriate for the individual convergence 
estimates Kest,i instead of a common residual distribution 
PA r( . s (Arcs). As Fig.plshows, the conditional distribution for 
an estimated convergence K cst — —0.05 is very narrow with 
a dispersion that is only 45% of the dispersion of the dis- 
tribution of the magnification over the full sky. A standard 
siren with a convergence estimate that low may thus provide 
substantially tighter constraints on cosmological parameters 
than a standard siren without a convergence estimate. On 
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Figure 6. The probability distribution of the residual magni- 
fication A rea = A — A est for sources at redshift z§ = 1.5. Com- 
pared are the distribution for A cs t = (solid line), the distri- 
bution PA roB (^ros) for the optimal estimate l |25| with K cs t recon- 
structed using the shear from an advanced survey with n ga i = 
lOOarcmin -2 and a filter scale 8 S = 25arcsec (dashed line), and 
the conditional distribution Px roa |K oat (-^rcsl^est) for K cs t = —0.05 
(dotted line) and re oa t = 0.1 (dash-dotted line). 



the other hand, a convergence estimate K cst =0.1 indicates a 
very broad magnification distribution, which yields a resid- 
ual dispersion that is twice the dispersion of the full-sky 
magnification distribution. Knowledge of such a high conver- 
gence estimate for a standard siren may help to avoid biases 
and underestimated errors on its distance due to lensing. 

The dispersion o"A rcs | reoat of the conditional probabil- 
ity distribution P\ TeB \ KeBt (A rea |>vest) as a function of the esti- 
mated convergence K cst from an advanced galaxy-shear sur- 
vey is shown in Fig. W\ This dispersion is the same for both 



the simple estimate ( |24[ ) and the optimal estimate ( 25 I (see 
AppendixlA} . The simple estimate, however, produces a con- 
ditional bias (A rcs ) A | , which is absent (by construction) 
for the optimal estimate. 

The small residual dispersion in regions with low values 
of the estimated convergence appears much more encourag- 
ing than the all-sky average residual dispersion (compare, 
e.g., the 45% residual for n cst — —0.05 to the all-sky aver- 
age residual of 78% a shear survey). The question is, how 
abundant are regions with small residual dispersion. Fig- 
ure [8] shows that the distribution of the conditional residual 
dispersion is very skewed. Values substantially below the 
all-sky average are common. For example, the median is at 
°"A r< , B |K CBt = 0.05 for an advanced shear survey, which is 10% 
below the all-sky average. 

As Fig.[9]illustrates, there is little benefit from the mag- 
nification correction for standard candles/sirens at redshifts 
zs much smaller than the median redshift z mc di an = 1-5 of 
the weak-lensing survey. For zs J; 0.5, a substantial part of 
the reconstructed convergence is due to matter structures 
at or above Zs- These matter structures shear most galaxy 
images used for the convergence reconstruction, but do not 
(de-)magnify the standard candle/siren, and thus merely act 
as noise in the magnification estimation. 

In contrast, standard candles/sirens at redshifts zs > 
^median do benefit from delensing. The residual dispersion 
a\ for sources at zs =3.1 can be reduced to 85% of the 
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Figure 7. Bias (A rcs ) 



Are 



and dispersion cr^ 



of the 



3 | Kest (Ares|«est) of the 
for sources at redshift 



conditional probability distribution p^ 
residual magnification A ros = A — A os t 
zs = 1.5 as a function of the estimated convergence K es t- Shown 
are the conditional bias and dispersion for the simple estimate 
( |24| with the convergence re cs t reconstructed using the shear from 
an advanced survey with n ga i = lOOarcmin -2 and a filter scale 
6 B = 25 arcsec. Note that the conditional bias of the optimal esti- 
mate ( |25| vanishes by construction, while its conditional disper- 
sion is equal to that of the simple estimate. 
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Figure 8. The probability distribution of the conditional resid- 
ual dispersion 0"a Ik * f° r sources at redshift z$ = 1-5 for an 
advanced survey with shear (dashed line). The vertical solid line 
indicates the value of the uncorrected magnification dispersion 
<j\, and the vertical dotted line marks the all-sky average <7\ 
of the residual dispersion. 
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Figure 9. The dispersion <TA rca of the residual magnification 
as a function of the source redshift zg. The dispersion without 
correction (Aest = 0, solid line) is compared to the residual dis- 
persion with convergence reconstructed from the shear (dashed 
line) and shear and flexion of an advanced weak-lensing survey 
with galaxy number density n ga i = lOOarcmin -2 and median 
redshift z mcdian = 1.5. 
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Figure 10. The probability distribution of the magnification- 
induced distance error 5Di um relative to the true distance D\ unl 
for sources at redshift zg = 1.5. Compared are the distribu- 
tion without correction (solid line), the all-sky distribution for 
the corrected distance using the optimal estimate ( |25| with K e st 
reconstructed using the shear from an advanced survey with 
n gal = lOOarcmin -2 and a filter scale a = 25arcsec (dashed 
line), and the distribution for sources in regions with estimated 
convergence K eB t = —0.05 (dotted line). 



dispersion ax for the uncorrected magnification, if only shear 
information is available. If flexion can be used in addition, 
the dispersion can be further reduced to 80% of the uncor- 
rected dispersion. 

Even for standard candles/sirens at redshift zs = 5.3, 
there is a noticeable reduction in the residual dispersion. The 
magnification correction based on shear yields a residual dis- 
persion o"A rcs = 0.9cta for the best filter scale 8 B = 25 arcsec. 
The best correction using shear and flexion reduces the resid- 
ual dispersion to 85% of the uncorrected dispersion. 

The distributions of the lensing-induced distance errors 
5-Dium are shown in Fig. [To] for sources at zs = 1.5. Before 
correction, only 50% of the standard candles/sirens have a 
relative distance error ^ 2%. After correction with shear 



Table 1. The fraction of sources at redshift zs with relative 
lensing-induced distance error e = |<5Di um /_Di um | below a given 
threshold before/after correction with shear data from an ad- 
vanced survey. 
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Figure 11. The dispersion &\ Ias of the residual logarithmic 
magnification A rcs = A — A cs t of sources at redshift zs = 1.5 
as a function of the filter scale B used to reconstruct the con- 
vergence Kest from a futuristic survey with n Ra ; = 500arcmin — 2 . 
Results are shown for the optimal estimate ( |2")| with convergence 
reconstructed from shear (dashed line) and from shear and flexion 
(dotted line). The solid horizontal line indicates the dispersion of 
Arcs = A without correction. 



data from an advanced survey, the fraction rises to 69%. In 
regions with K eBt = —0.05, 90% of the sources have a lensing- 
induced distance error below 2%. The fraction of sources 
with absolute lensing-induced distance error below a given 
threshold before and after correction is listed in Table [l] for 
several source redshifts. 



4.3 Futuristic reconstruction 

Improvements in observation instruments and techniques 
might one day permit weak-lensing surveys with very high 
galaxy number densities. As a futuristic scenario, we con- 
sider a very deep high-resolution survey that provides shear 
and flexion estimates from galaxies with a density n ga i = 



500 arcmin and median redshift Zmedia 



1. 



As Fig. [TT] shows, the high number density allows one 
to substantially reduce the residual dispersion (TA ro3 for stan- 
dard candles/sirens at zs = 1.5 at small smoothing scales 8 S . 
If only shear data are used, the residual dispersion <j\ ios can 
be reduced by delensing to 68% of the uncorrected disper- 
sion a\ at 8 S = 10 arcsec. If shear and flexion data are used, 
the residual dispersion is 0.04 at 9 S = 3 arcsec, which is only 
54% of the uncorrected dispersion. To reach a similar statis- 
tical uncertainty due to gravitational magnification without 
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Figure 12. Dispersion cr\ It , s i Kost of the conditional probabil- 
ity distribution P\ IOB \ K „t (A re s|K e st) °f * ne residual magnification 
Arcs = A — Aest for sources at redshift zg = 1.5 using the optimal 
estimate \25\ and convergence K cs t reconstructed using the shear 
and flexion from a futuristic survey with n ga i = 500 arcmin -2 and 
a filter scale 8 S = 3 arcsec (dashed line) . The solid line indicates 
the dispersion of the uncorrected magnification. 



0.15 



l 0.10 



0.05 



0.00 




zs 



Figure 13. The dispersion o"A rcB of the residual magnification 
as a function of the source redshift zg. The dispersion without 
correction (A cs t = 0, solid line) is compared to the residual dis- 
persion with convergence reconstructed from the shear (dashed 
line) and shear and flexion of a futuristic weak-lensing survey 
with galaxy number density n ga i = 500 arcmin -2 and median 
redshift z mcdian = 1.8. 



delensing, one would need at least a four times higher num- 
ber of similar standard candles/sirens. 

For low values of the estimated convergence n cst , the 
residual dispersion is even lower, as can be seen in Fig. |12| 
For example, the residual dispersion o"A rcs |fc oat f° r K cst = 
—0.05 is only 25% of the unconditional uncorrected disper- 
sion a\. 

Figure [13] illustrates that a futuristic survey may also 
help to substantially reduce the magnifications scatter for 
high- redshift standard candles/sirens. For example, the 
residual dispersion for sources at redshift 2s = 2.1 can be 
reduced to 50% of the uncorrected dispersion with a futur- 
istic survey measuring shear and flexion. This is even better 
than for the sources at redshift zs = 1.5. For zg = 5.3, the 
residual dispersion can be reduced to 67% of the uncorrected 
dispersion. 
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Figure 14. The probability distribution of the magnification- 
induced distance error SD\ unl (relative to the true distance Di um ) 
for sources at redshift zg = 1.5. Compared are the distribution 
without correction (solid line), the all-sky distribution for the 
corrected distance using the optimal estimate i ]2"i| with re os t re- 
constructed using the shear and flexion from a futuristic survey 
with n ga ] = 500 arcmin" 2 and a filter scale 8 S = 3 arcsec (dashed 
line), and the distribution for sources in regions with estimated 
convergence K cs t = —0.05. 



Table 2. The fraction of sources at redshift zg with relative 
lensing-induced distance error e = |<5-Di um /-Dluml below a given 
threshold before/after correction with shear and flexion data from 
a futuristic survey. 







e < 0.01 


e s£ 0.02 


e < 0.05 


e s£ 0.10 
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= 1.0 


0.41/0.64 


0.80/0.88 


0.96/0.98 
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ZS 
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ZS 


= 4.2 


0.11/0.21 


0.22/0.41 


0.54/0.81 


0.88/0.97 


zs 
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0.49/0.73 


0.83/0.94 



Fig. 1 14| shows for sources at zs = 1.5, how well a futur- 
istic survey with shear and flexion data can help to mitigate 
the lensing effects on distance errors. While the fraction of 
uncorrected standard candles/sirens with relative distance 
error ^ 2% is only 50%, the fraction of corrected standard 
candles/sirens is 82%. If only regions with K ost = —0.05 are 
considered, 97% of the sources have residual distance er- 
rors ^ 2%. Results for several distance error thresholds and 
source redshifts are listed in Table [2] 



4.4 Improving the reconstruction with Wiener 
filters 

Using a Gaussian filter to smooth the reconstructed conver- 



gence [cf. Eq. ( 23 1] might not be the best way to exploit the 
spatial correlation in the convergence field. With a few as- 
sumptions, choosing the optimal filter for the reconstructed 
convergence becomes a text-book problem on Wiener filters 
( |Wiener|1949| . 

Consider the unfiltered convergence re raw reconstructed 
from galaxy shear or shear and flexion as a raw estimate for 
the convergence k(zs) to a standard candle/siren at redshift 
zs- Let -P Kraw (^) denote the power spectrum of K raw , and 
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Figure 15. Wiener niters W(£) in Fourier space for stan- 
dard candles/sirens at redshift 2g = 1.5 and convergence re- 
constructed from the shear (solid line) or the shear and flex- 
ion (dashed line) from an advanced survey with galaxy density 
n gal = lOOarcmin - 2 and median redshift z mc dian = 1-5 (a) and 

^median t.o 



a futuristic survey with n 

(b). 



gal 



500 arcmin and 



-P«;(z s ),«; raw W denote the cross power spectrum of k(zs) and 
Kraw Then 



W(£) = 



«*S),' 



>(*) 



AQ 



(27) 



provides a Wiener filter in Fourier space, which minimises 
the mean square difference between the 'true' convergence 
k(zs) and estimated convergence K eat (0), given in Fourier 
space by 



ft-tCe) = w(\£\)k^ w (£). 



(28) 



The Wiener filter ( 27 1 can be constructed by measuring 



the power spectra in the simulated fields and approximating 
their ratios with a suitable function (e.g. a spline). The niters 
for standard candles/sirens at redshift zs = 1.5 are shown 
in Fig. [15] for the advanced and the futuristic survey. 

Using Wiener filters for the estimated convergence im- 
proves the dispersion in the residual magnification only 
marginally compared to the best Gaussian smoothing. For 
an advanced survey with n ga i = lOOarcmin -2 with shear 
data, the Wiener filter yields a residual dispersion crA rcB = 




-0.2 0.0 0.2 

nfe = 3.1) 



Figure 16. The joint distribution of the shear component 71 
at redshift z = 1.8 and redshift zq = 3.1 at the same sky po- 
sition (lighter areas correspond to higher probability densities 
on a logarithmic scale). The dashed line marks the mean shear 
(7l( 2 )) 71 (z)| 71 (z s ) as a function of 7i(zs)- Note that the joint 
distribution for the second shear components 72(2) and 72 (zs) is 
the same as for the first shear components. 



78% of the uncorrected dispersion <j\, whereas the Gaus- 
sian filter with 9 E = 25 arcsec yields a residual dispersion 
°"a ios = 79%. For a futuristic survey with shear and flexion, 
the residual dispersion is 53% of uncorrected dispersion for 
the Wiener filter, which is to be compared to 54% for the 
best Gaussian filter. 



4.5 Improving the reconstruction with individual 
galaxy redshifts 

The image distortion for a galaxy at redshift z, and hence 
the information in the galaxy image about the magnifica- 
tion at the redshift zs of a standard candle/siren is not uni- 
form in redshift z. In a full tomographic analysis, the varia- 
tion of the lensing signal with redshift is used to perform a 
three-dimensional reconstruction of the matter structures in 
the light cone (e.g. Massey et al.| |2007). The reconstructed 
three-dimensional matter distribution can then be used to 
estimate the magnification for a given sky position and red- 
shift. Here, we take a simpler approach, which does not try 
to locate the actual matter structures causing the lensing 
in redshift, but only exploits the statistical relation between 
lensing quantities at different redshifts. 

We start from the same assumptions used for the simu- 
lated the weak-lensing surveys (see Section [3. 2\ : The image 
shape of a galaxy at sky position # ga i and redshift z ga i pro- 
vides noisy but unbiased estimates 7 ga i , J-" ga i , and C/gai of the 
shear 7(0 ga i,z g ai) and flexion T{0 ga \, z ga i) and <5(0 ga i, z ga i) 
to the galaxy redshift z ga i. The error in these estimates is 
simply the shape noise with vanishing mean. Considered as 
estimates for the shear and flexion 7(0 ga i,zs), F{Q sa .\,zs), 
and (/(0 g ai, zs) at the redshift zs of a standard candle/siren, 
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H 0.8 ; 




Figure 17. The ratio r 7 (z,zs) between the mean shear at red- 
shift z and the shear at redshift zg = 3.1. The ratio has been 
calculated by fitting the relation | |29a| l to the shear data from the 
ray-tracing simulations. 



7gai, ■T'gai, and <5 ga i acquire an additional error, namely the 
deviation between the shear and flexion at different red- 
shifts. In the following discussion, this deviation is split into 
a systematic (i.e. zero-scatter) and a statistical (i.e. zero- 
mean) component. 

The first step towards an improved reconstruction aims 
at eliminating the systematic deviation. As Fig. [16] illus- 
trates, the shear values for two different redshifts along the 
same line of sight are strongly correlated, but on average, 
the shear to the lower redshift is smaller in magnitude. Fig- 

of the shear 



16 



also shows that the mean (■y(z)) , ,, , 
y(Z) to redshift z for given shear •y(zs) to redshift zs along 
the same line of sight is well approximated by the linear 
relation 

<7(z)} 7(z) , 7(zs) = r 1 {z, zsh(zs) (29a) 

with a redshift-dependent proportionality factor r 7 (z, zs) 



(see Fig. 17 1. Moreover 



(^( z )) >( Z )|.F( ZS ) =rr(z,zs)J 7 (zs) and 
(0( 2 )>eO)ie(zs) = r e(z,z s )G(z s ) 



(29b) 
(29c) 



with factors rj?(z, zs) and rg(z, Zs) very similar to r 1 (z, zs). 
Thus, one can construct unbiased estimates 7 ga i(«s), 
•^gai(zs), and Q sa ,i(zs) of the shear and flexion at redshift 
zs by: 



7gai(zs) = [r 7 (2 gB i,2 S )] _ 7gal, 
■Fgai^s) = [rr(z sa i, zs)] T gaU and 

S K al(Zs) = Wg(z e& l,Zs)] C/gal. 



(30a) 
(30b) 
(30c) 



In a second step, the statistical weights of the galax- 
ies in the reconstruction are extended to take into account 
the redshift dependence of the error in the scaled estimates 
(301. The error in these estimates has two components: the 
shape-noise and the scatter around the mean relation be- 
tween the lensing quantities at different redshifts. The shape 
noise in the scaled estimate is simply the shape noise in 
the unsealed estimate (de-) amplified by the proportional- 
ity factor r.y(z,zs), rj:{z,zs), or rg(z, zs), and has no spa- 
tial correlation prior to smoothing. The scatter around the 
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Figure 18. The dispersion crx Iog of the residual magnification as 
a function of the source redshift zg . The dispersion without cor- 
rection (A cs t = 0, solid line) is compared to the residual disper- 
sion with convergence reconstructed from the shear and flexion 
of a futuristic weak-lensing survey with galaxy number density 
n gal = 500arcmin — 2 and median redshift z mcc iian = 1-8, when 
ignoring (dashed line) or using the individual galaxy redshifts 
(dotted line). 



mean relation between shear and flexion at different red- 
shifts has a quite different redshift dependence (e.g. it van- 
ishes when Zgai = zs), and has a finite spatial correlation 
inherited from the intrinsic spatial correlation of the shear 
and flexion at a single redshift. As a consequence, the opti- 
mal redshift weights depend on the galaxy number density 
and the spatial smoothing. However, an examination of the 
optimal weights (which can be derived in a lengthy calcula- 
tion) reveals that the contribution from the scatter around 
the mean relation is much smaller than the shape-noise con- 
tribution for the smoothing scales and number densities of 
interest (at least for z < zs). Thus one can approximate the 
redshift weights by: 



W z ,~,(z,Zs) = r- 1 (z, z s ) , 
W z ,f(z, z s ) = r F (z, zs) , and 
w z ,g(z,zs) — rg(z,z S ) . 



(31a) 



Optimal use of the scaled estimates ( 30 1 and redshift- 
dependent weights ( |31| l requires that the information about 
the error of the scaled estimates is still available for the 
spatial smoothing of the resulting fields. Thus, the recon- 
struction algorithm discussed in Section |3.3|is modified by 



merging the spatial-filtering step with the step where the 
galaxy estimates are projected onto a regular mesh. Since 
the smoothing is now computed in real space, the algorithm 
is much slower than the one described Section |3.3| in par- 
ticular for larger smoothing scales. 

For standard candles/sirens at redshift zs = 1.5, there 
is no improvement in the reconstruction from an advanced 



or futuristic survey (shown in Fig. 18 1. This is different for 
higher redshifts. For zs = 3.1, the reconstruction from the 
shear and flexion of a futuristic survey changes from 57% 
to 52% of the uncorrected dispersion if redshift weights are 
used. For zs = 5.3, redshift weighting improves the residual 
dispersion from 67% to 62%. 
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5 SUMMARY AND DISCUSSION 

Unless corrected for, gravitational lensing induces significant 
errors in the measured distances of high-redshift standard 
candles and standard sirens. In this work, we have used nu- 
merical simulations to investigate how much these distance 
errors can be reduced with weak-lensing reconstruction. The 
method considered comprises (i) reconstructing the conver- 
gence towards a standard candle/siren from the shear and 
flexion measured in a galaxy lensing survey; (ii) estimat- 
ing the magnification from the reconstructed convergence; 
and (iii) correcting the observed signal of the standard can- 
dle/siren using the inferred magnification. 

By measuring the relation between the reconstructed 
convergence and the magnification, we have constructed an 
optimised magnification estimate that is unbiased and min- 
imises the residual magnification errors. Furthermore, we 
have studied the optimal smoothing scale for the weak- 
lensing reconstruction under various survey conditions. 

For an advanced shear survey with n ga i = 100 arcmin -2 , 
we find that the lensing-induced distance errors for standard 
candles/sirens at redshifts zs ~ 1.5 are reduced on aver- 
age by 20%. This confirms earlier findings by |Dalal et al.| 
(2003) that reconstructions based on such surveys do not 
significantly reduce the lensing errors on average. However, 
in regions with low estimated convergence n CBt , the residual 
magnification uncertainty is substantially smaller than in re- 
gions of high convergence (e.g. only 40% of the uncorrected 
error for K est = —0.05 compared to 200% for « cst = 0.1). 
Thus a weak-lensing survey is useful to identify standard 
candles/sirens in regions with low (high) convergence and 
small (large) magnification errors, which can then be given 
larger (smaller) weights in a statistical analysis. 



As already pointed out by Shapiro et al. ( 2010 1, the cor- 



rection from weak-lensing reconstruction can be improved 
considerably by including flexion measurements and greatly 
increasing the galaxy number density. For example, a fu- 
turistic shear and flexion survey with n ga i = 500 arcmin -2 
yields error reductions of 50% for standard candles/sirens at 
zs ~ 1.5, and 35% for sources at zs ~ 5. In low-convergence 
regions, the errors can be reduced by up to 75%. For exam- 
ple, the residual distance errors are below 2% for over 97% 
of the standard candles/sirens in regions with K cs t = —0.05. 
Such lensing-corrected standard candles/sirens would con- 
stitute very competitive distance indicators (but one should 
keep in mind that these numbers are based on very opti- 
mistic assumptions about the properties of future surveys). 

Our simulation provides detailed statistical information 
about the magnification-induced distance errors before and 
after correction, in particular their probability distributions 
as a function of the estimated convergence. The distributions 
are valuable for Bayesian parameter estimation and model 
selection. In a forthcoming paper (Gair, King, & Hilbert, in 
prep.), we will discuss the resulting implications for the use 
of high-redshift standard sirens as cosmological probes. 

In this paper, we have considered several ways to im- 
prove the accuracy of the magnification correction scheme. 
Wiener filters constructed from the measured (cross-)power 
spectra of the true and reconstructed convergence are ex- 
pected to provide the optimal filtering of the noisy conver- 
gence maps. However, in our simulations, they perform only 



marginally better than Gaussian filters with suitable filter 
scale. 

Furthermore, we have investigated a simple method em- 
ploying redshift-dependent weights for the galaxies to im- 
prove the weak-lensing reconstruction. While standard can- 
dles/sirens at redshifts zs ^ 2 show no benefit from this 
method, magnification errors for higher redshifts are reduced 
by an additional 5%. Further improvement might come from 
a tomographic reconstruction (Simon et al. 2009). 

In contrast to correction schemes based on modelling 
the foreground matter structures from the observed light 
(e.g 



Gunnarsson et al. 20061, the method discussed here 



has the advantage that it does not require any assumptions 
about the relation between visible and dark matter. On the 
other hand, the light distribution contains valuable informa- 
tion about the matter distribution and the magnification. A 
better magnification correction might be obtained by infer- 
ring the foreground matter distribution from both the ob- 
served properties of galaxies along the line of sight and the 
shear and flexion information from distant galaxies. Addi- 
tional information about the matter distribution at redshifts 
beyond those probed by conventional galaxy weak lensing 
and flexion could be gathered from observations of lensing 



of the cosmic microwave background (Lewis and Challinor 
|2006j ) or high-redshift 21-cm radiation ( Hilbert et al.|2007a[ ) 

In our simulations, we have fully taken into account 
many aspects of weak-lensing galaxy surveys that impact 
the accuracy of the correction. These include, e.g., the non- 
Gaussian nature of the convergence field, the randomness 
and discreteness of the galaxy redshifts and image positions, 
and the complex relation between the reconstructed conver- 
gence and the magnification. Thus, our studies provide more 
realistic predictions about the accuracy of the correction 
than earlier studies relying on convergence power spectra. 

However, our approach neglects several complications 
that affect the correction. Future studies should include the 
effects of observing reduced shear and flexion. The impact of 
source clustering and intrinsic alignment on the weak-lensing 
reconstruction must be investigated. 

Our simulations do not include the effects of lensing 
by structures on scales below the resolution of the MS (i.e. 
dark-matter structures on scales < 5/i -1 kpc or luminous 
structures with masses < 10 9 /i -1 M©). Such structures could 
affect the magnification of SNe, GRBs, and SMBs, but it 
would be difficult to recover them in a weak-lensing mass re- 
construction. The resulting increase in the uncorrected and 
residual magnification error should be investigated, possi- 
bly with higher-resolution simulations. One should also con- 
sider using the observed scatter in samples of standard can- 
dles/sirens to constrain the small-scale matter fluctuations 
missed by the simulations. 

Finally, the cosmology dependence of the correction 
should be explored. For example, a higher/lower normali- 
sation as increases/decreases both the uncorrected and the 
residual magnification dispersion. Moreover, errors in the 
assumed cosmology lead to a suboptimal magnification es- 
timator and thus increase the residual magnification error 
and might also introduce a bias. The resulting degradation 
should be quantified in future studies. One should also in- 
vestigate how well one could detect a possible bias in the 
magnification estimates by looking at differences between 
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the expected and observed distributions of the convergence 
and magnification estimates. 
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APPENDIX A: MORE ABOUT 
MAGNIFICATION ESTIMATES FROM NOISY 
CONVERGENCE MAPS 

Define for any two random variables x and y with joint prob- 
ability density function (pdf) Px, y (x,y) and conditional pdf 
p x \ y (x\y), and any function f(x,y) as short-hand notation 



residual equals the conditional variance a x 



of the true 



dxdyp x , y (x,y)f(x,y) 



if) = </>„,„ : 

for the unconditional mean of /, and 

(f)\ v (y) = (f} x \ y (y) = dx Px\y( x \y)f( x ,y) 



(Al) 



(A2) 



for the conditional mean of / for given y. Furthermore, de- 
fine 



(A3) 



°7 = dxd yp*,y( x ,y)[f( x ,y) - (/>]' 

for the full variance of /, and 

°}\v{v) = dx Px\y( x \y)[f( x ,y) - (f) x \ 



(A4) 



for the variance of / restricted to a fixed y. 

Consider an ensemble of lines of sight (l.o.s.), where 
each l.o.s. is characterised by its (logarithmic) magnification 
A and its estimated convergence K cst . Assume that the joint 
pdf p\, Kast (A, /test) of A and K ost is known for the ensemble, 
which implies that the marginal and conditional distribu- 
tions of A and K cst are known, too. Furthermore, assume the 
convergence K cst known for each l.o.s. 

Consider a magnification estimate A os t = A cst (K e st) that 
is based on the estimated convergence K e st- Consider the 
(unknown) residual magnification 



t(^est) 



(A5) 



obtained after correcting the (unknown) magnification A of a 
l.o.s. by the magnification estimate A cs t(Kcst) obtained from 
its (known) convergence estimate K cat . The conditional bias 
(Ares) \ i , i-e. the mean residual for a fixed estimated con- 
vergence K os t, is then given by: 



(^ rcs )A| Kcat ( K <=st) — W\\ Kost («est) - A ost (K cst ). 



(A6) 



Moreover, the conditional distribution PA rcB |K oat (Arcsl^est) of 
the residual A res for a given convergence estimate K est is then 
simply a shifted version of the conditional pdf PA|K C3t (A|«est) 
of the true magnification A: 



PA rca |fc oat (Arcs|K cst ) — P\\ KcBt (Arcs + Aest (fvest ) | Ke 

This implies that the conditional variance <r\ 
© 0000 RAS, MNRAS 000, 000-000 



sliest 



(A7) 
of the 



magnification for any magnification estimate A C st(«est). 

The problem of finding the optimal magnification es- 
timate A st(Kcst) can be approached with basic calculus 
of variations. First, one has to define 'optimal', e.g. by 
demanding that an optimal estimate (i) is unbiased, i.e. 



(Aest); 

an«El 



w 



A|k c 



and (ii) minimises the residual vari- 



2 
°"A r , 



Then, 



= 



c) 



d\ 



A\AK cst P\^ ast {\, K cst )[\ — Aest(Kest)]" 
dKostP Kost ( K cst) 
X /dAp A | Kcat (A|K C st)[A - A oat (K est )] . 

- d\p x \ Kcst (A|K cst ) [A — Aest(«est)] 



(A8) 



'dA PA | Kcst (A| Kost )[A-A C st(«;est)] (A9) 

= (^)A| Kcst ( K <=st) — A cst (K cst ) 

for the optimal estimate. Thus, the optimal estimate reads: 

Aest(Kest) = W \\ Kcst («est). (A10) 
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